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This article presents a newly developed modification of the dandelion 
optimizer (DO). The proposed method is a chaotic algorithmic integrity and 
modification of the original dandelion optimizer. Dandelion is one of the 
plants that rely on wind for seed propagation. This article presents the tuning 
of the power system stabilizer with the method proposed in a case study of a 


single machine system. The validation of the proposed method uses the 


benchmark function and performance on a single engine system against 
Keywords: transient response. The method used as a comparison in this article is the 
whale optimization algorithm (WOA), grasshopper optimization algorithm 
(GOA) and the original dandelion optimizer (DO). The simulation results 
show that the proposed method, which is a modified dandelion optimizer, 
provides promising performance. 
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1. INTRODUCTION 

The use of renewable energy and distributed generation integrated with various power semiconductor 
devices affects the stability of the power system [1]—[4]. The robustness of the interconnection of the electric 
power system will decrease because it is influenced by low-frequency oscillations that are not optimally 
submerged [5], [6]. This will cause electricity distribution disturbances which will indirectly burden the system 
technically and economically. So, the focus on attenuation of low-frequency oscillations becomes an important 
key in power systems. 

A popular device in dealing with low-frequency oscillations is the power system stabilizer (PSS). This 
device is used in a variety of system conditions. Conventional PSS is modeled linearly as the basis of its design. 
This model considers the optimal working point as the basis for tuning PSS parameters. An electric power 
system that has non-linear and variable characteristics over a wide range. This makes conventional PSS not 
optimal. The development of technology and computational algorithms that are developing rapidly. This affects 
and penetrates into the power system stabilizer. Several studies have presented findings regarding this matter. 
In decades, computational approaches have been presented in tuning PSS parameters such as atomic search 
optimization (ASO) [7], moth search algorithm (MSO) [8], whale optimization algorithm (WOA) [9]-[12], 
Henry gas solubility optimization (HSO) [13], particle swarm optimization (PSO) [14]-[17], grey wolf 
optimization (GWO) [18]-[21] and JAYA algorithm (JA) [22]—[24]. 

This article presents the power system stabilizer tuning method using the modified dandelion 
optimizer (DO) method. Dandelion is one of the plants that rely on wind for seed propagation [25]. DO 
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modification method is to integrate the chaotic algorithm and change one of the DO components. This is as a 

goal to increase the ability of DO. The contributions of this research are, 

— Modify the dandelion optimizer method called CDO to get a new balance point of exploration and 
exploitation. 

— Application of the CDO method to tune PSS parameters by measuring performance by comparing with 
conventional methods (PSS-Conv), PSS based on whale optimization algorithm (PSS-WOA), PSS based 
on grasshopper optimization algorithm (PSS-GOA) and PSS based on dandelion optimizer original 
(PSS-DO). 

This article is structured: The second part describes the dandelion optimizer method, the dandelion 
modification method and the power system stabilizer. The third section presents the results and analysis. The 
last section provides conclusions from the research. 


2. METHOD 
2.1. Dandelion optimizer (DO) 

Dandelion optimizer (DO) is a method adopted from the movement of plant seeds. Dandelion is one 
of the plants that rely on wind for seed propagation. Two important factors that affect the spread of dandelion 
seeds are wind speed and weather. The falling distance of dandelion seeds is affected by wind speed. 
Meanwhile, weather affects the ability of seedlings to grow near or far. DO can be modeled mathematically in 
3 parts, namely the ascending section, descending section, and landing section. DO is like a population-based 
algorithm that assumes every dandelion seed is a candidate solution. 


xt wee xe 
A=|i =: (1) 
ey eee a 
Xi =Xmin + rand (Xmax = Xmin) (2) 
fp = min (f(x;)) (3) 
Xe = x(find (fy == f(%))) (4) 


Where population size is symbolized by p. Variable dimensions are represented by d. rand denotes a random 
[0,1]. 


2.1.1. Ascending section 

In the rising stage, sunny and windy weather brings dandelion seeds up. On the other hand, there is 
no wind over the seed when it rains. Local model search occurs in this section. Flying dandelion seeds are 
affected by wind speed and humidity. Dandelion seeds have the characteristic of being able to fly far 
considering the height. In this section the weather is modeled in two namely. 

Condition 1: Sunny day conditions, lognormal distribution applied as wind speed. In this session, DO 
conducts exploration. The wind causes dandelion seeds to move randomly to various locations. The wind speed 
determines the height of the seed. Session 1 can be modeled, 


Xpay =X +aXv, xX vy XINY X (x, — x) (5) 
x, =rand (1,dim ) X rand (X%max — Xmin) + Xmin (6) 
=< a 2 
Iny = yuan Plage On) ly = 0 (7) 
0 y<0 
a = rand ()x (St? -=t+ 1) (8) 
Vy =1T X cosé (9) 
vy =r xsind (10) 
r= = (11) 
6 =(2 Xrand()-1)xza (12) 
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Where the position of the dandelion seed during iteration is symbolized x,. Randomly selected 
position in the search space during iteration is symbolized x,. InY is a lognormal distribution. The adaptive 
parameter used to adjust the search step length is a. The coefficient of the dandelion rising passage as a result 
of the action of the separate eddies symbolized v, and vy. 

Condition 2: that is on a rainy day. Dandelion seeds have problems growing. In this condition, local 
exploitation is carried out. The mathematical equation of condition 2 is, 


Xpa1 =X XK (13) 
| 2__ 2 a 

d= Flore T2-2T+1 a T2-2T+1 (14) 

k =1-randQ xq (15) 


os = aan ae x In¥ x (%; — X) randn < 1.5 (16) 


ua Xk else 


where k is applied to maintain the local seeking domain of an agent, randn represented the random value that 
obeys the basic normal distribution. 


2.1.2. Descending section 

At this stage, the exploration phase is emphasized. The movement of dandelion seeds will decrease 
for sure after experiencing a peak at a certain value. The average information after the ascending stage is used 
to reflect the stability of the parental offspring. This is to provide support for the improvement of the overall 
population. The mathematical modeling of this stage is, 


Xta1 = Xp — AX Be X (Xmean_t — & X Be X Xt) (17) 
1 
Xmean.t = sop List Xj (18) 


where f£, Points Brownian action. It is a random value from the standard normal distribution. 


2.1.3. Landing section 

Exploitation phase occurs in this section. The landing place of the dandelion seeds is chosen at 
random. The approximate position of the most viable dandelion seed was used as the optimal solution. Elite 
information is currently exploited in the local environment to obtain global optimum accuracy. This behavior 
can be modeled, 


Xt=1 = Xetite + levy(A) X @ X (Xelite — Xt X 6) (19) 
levy(A) = 5°? (20) 
alg 

r(+p)xsin@4) 

OS pate (21) 
r@t®)xax202 > 

= 
6=3 (22) 


where Xerite reflects the best position of the agent in every iteration. levy(A) deputizes the function of Levy 
flight. 


2.2. A novel modified dandelion optimizer (CDO) 

The proposed algorithm is an integration of the chaotic algorithm (CO) and the modified Dandelion 
Optimizer. This method is proposed to improve the DO optimization algorithm. CO applies chaotic variables 
instead of random variables. Chaos has non-reinforcement and ergodistic characteristics. In addition, the search 
system has a higher speed compared to search methods that are stochastic or rely on probability [26]. This 
study uses 1-D non-reversed maps, namely logitstics as a chaotic set algorithm. Modification is used to 
accelerate the level of the convergence curve to reach the optimal point. 


ylog(is1) = aX ylogyy(1 = ylogiiy) (22) 
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In (22) is used to replace the variable rand () in (12). So, (12) becomes, 

0 = (2X yloggi41) —1) X 0 (23) 
the second step is to change (19) to (24). This modification aims to sharpen the convergence curve obtained. 

Xt=1 = Xelite + levy(A) X @ X (Xetite — Xt X 5) — Xelite X rand() (24) 


2.3. Power system stabilizer 

The damping torque on the engine rotor will be regulated by PSS with the aim of producing a 
compensation between the electric torque and the excitation input [27]. PSS will output a value proportional to 
the rotor speed. This maintains the stability of the electrical system. Figure | is a lead-lag PSS schematic. 


2.4. Design of controllers 

CDO is used to find the optimal point of attenuation by adjusting the PSS parameter. So that the 
requirements of the transient response criteria can be increased in the closed loop response. An illustration of 
PSS tuning with CDO can be seen in Figure 2 on Appendix. The initial step starts with making the required 
parameters. The results obtained are random parameters which are always corrected during the iteration process. 
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Figure 1. PSS lead-lag type 
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3. RESULTS AND DISCUSSION 
3.1. Convergence curve profile 

The performance of the CDO Algorithm is measured using the classical and CEC 2019 benchmark 
function. In this article, 12 classical benchmark functions are used which can be seen in detail in Table 1. The 
results of the classical benchmark function on CDO compared to the GOA, WOA, and DO methods can be 
seen in Figure 3. The Unimodal functions can be seen in Figure 3(a) to Figure 3(d), Multimodal can be seen in 
Figure 3(e) to Figure 3(h) and multimodal with fixed dimensions is presented in Figure 3(i) to Figure 3(1). The 
result of classical benchmark function can be seen in Table 2. This article also tests using CEC2019. The test 
results with CEC 2019 can be seen in Figure 3(m) to Figure 3(n). The Detail of CEC2019 Benchmark Function 
can be seen Table 3. From the results of trials using several functions from the CEC 2019 Benchmark function, 
the CDO method has a better average convergence value than the DO, WOA and GOA methods. The results 
of CEC 2019 Benchmark Function can ben seen in Table 4. 


Table 1. The classical benchmark functions 


ID Test Function Range Type 
n Unimodal 
ie F,(x) = yx [-100,100]" 
i=1 
ue Ua Unimodal 
F, Fe) => al +] [lx [-10,10)" 
i=1 i=1 
ae Unimodal 
F; F3(x) =O [-100,100]” 
i=1 j=1 
F, F,(x) = max {|X;|,1 <i<n} [—100,100]" Unimodal 
us Multimodal 
Fy Reo >. —X, sin(/ XI) [-500,500]" 
i=1 
n Multimodal 
Fy F(x) = Ye — 10 cos(2nX;) + 10] [-5.12,5.12]” 
i=1 
Multimodal 
Fi F,9(x) = —20 exp(—0.2 cos (27X;)) +20 +e [-32,32]” 
Multimodal 
Fis F(x) = wan * ef fads a [-20,20]" 
Fyg(x) = [1 + (4, +x. + 1)?19 - 14x, + 3x,2 — 14x, + 6x,x, + 3x,7)] x [30 Multimodal with 
Fy, + (2x, [-2,2]? fixed dimensions 


— 3x2)? x 8 — 32x, + 12x,? + 48x, — 36x,x2 + 27x,7)] 
Multimodal with 


Fyo Fig(x) = >) C; exp( — x aij (Xi — pij)”) [0,1]8 fixed dimensions 
Multimodal with 
Foo Fy9(x) = >) C; exp( — x aij (Xi — pij)*) [0,1]° fixed dimensions 


Multimodal with 
Foy Fy, (x) = Sia —a;)(X —a,)7 + CJ [0,10]* fixed dimensions 
j=1 


3.2. Transient response validation 

In this article, small signal stability analysis is applied to the Heffron-Philips model. The non-linear 
set of PSS parameters is solved by CDO. The parameters obtained are used to reduce wave oscillations as best 
as possible. This study uses a comparison method, namely conventional PSS, WOA, GOA and DO as a 
validation of the performance of the CDO. The PSS parameters that have been obtained are tested with the 
system under 100% loading conditions. Performance validation is carried out by comparing the CDO method 
with other algorithms and can be seen in Figure 4. The analysis of the transient response is detailed in Table 5. 
The simulation results using the MATLAB/simulink application with a laptop that has an Intel I5-5200 2.19 
GHz processor specification and 8 GB RAM memory where the PSS-CDO method is able to reduce overshoot 
from a speed of 96.48% from PSS-Conv. It can be seen in Figure 4(a). Meanwhile, the undershoot of the PSS- 
Conv method can be reduced by 75.99%, by implementing PSS-CDO. The Rotor angel comparison results can 
be seen in Figure 4(b). 
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Table 2. Results of the benchmark function 


ID CDO DO GOA WOA 
Fl Best 1.23E-11 642.4311 530.3973 0.001318 
Worst 3.25E-10 2404.2 1185.804 0.3713 
Mean 1.29E-09 = 4219.987 2217.168 3.4714 
Std Deviation 3.24E-10 982.7503 409.3191 0.67814 
F2 Best 3.16E-07 8.6321 8.567 0.005495 
Worst 4.75E-06 23.2267 27.9892 0.084016 
Mean 1.16E-05 76.2167 85.1458 0.27367 
Std Deviation 2.61E-06 13.9029 21.4767 0.074365 
F3 Best 1.95E-10 — 8200.342 1427.746 49283.88 
Worst 5.30E-09 16207.47 5290.209 100263 
Mean 2.52E-08 — 23853.42 15021.57 189167.9 
Std Deviation 7.40E-09  4707.366 3251.914 31384.76 
F4 Best 4.02E-06 35.8668 9.596 16.9667 
Worst 1.31E-05 49.8996 15.224 66.0116 
Mean 2.73E-05 69.733 23.7969 89.2632 
Std Deviation 5.96E-06 9.14 3.7136 20.5823 
F8 Best -3437.92  -8086.5344  -7407.9610 -11917.8502 
Worst -2491.49  -6574.3898 — -6332.4420 -9133.0283 
Mean -1578.35  -5285.6636 = -4765.8226 -7066.2870 
Std Deviation 469.98 770.5834 653.0869 1191.3498 
F9 Best 1.24E-11 74.7104 111.4599 0.0042695 
Worst 1.09E-10 132.3206 180.5123 54.3783 
Mean 3.59E-10 202.8456 243.5834 269.2653 
Std Deviation 8.60E-11 33.0798 29.7397 84.7164 
F10 Best 1.09E-06 8.3848 6.5252 0.010797 
Worst 3.41E-06 12.2815 8.9729 0.11812 
Mean 9.98E-06 16.0868 11.8267 0.43745 
Std Deviation 2.27E-06 2.5945 1.5671 0.12948 
Fil Best 8.27E-12 7.2046 4.1278 0.011716 
Worst 4.84E-10 21.5406 10.2899 0.28451 
Mean 2.08E-09 47.1931 19.3694 1.0244 
Std Deviation 5.79E-10 9.2638 3.5281 0.3107 
F18 Best 3.00 3.0000 3.0000 3.0000 
Worst 13.69 11.6400 3.0000 7.3665, 
Mean 92.20 84.0001 3.0000 30.6394 
Std Deviation 29.54 23.0161 0.0000 10.2109 
F19 Best -3.86 -3.8628 -3.8628 -3.8623 
Worst -3.76 -3.8627 -3.8037 -3.8290 
Mean -3.46 -3.8619 -2.9627 -3.6659 
Std Deviation 0.11 0.0002 0.1861 0.0544 
F21 Best -3.55 -10.1532 -10.1532 -9.8812 
Worst -1.97 -5.3342 -5.0476 -6.4502 
Mean -0.60 -2.6303 -2.6305 -2.5763 
Std Deviation 0.89 3.2221 3.0965 2.4101 


Table 3. CEC 2019 benchmark function [28] 


ID Function Range D 
Cec01 Storn ’s Chebyshev Polynomial Fitting [-8192, 8192] 9 
Cec02 Inverse Hilbert Matrix [ -16,384, 16,384] 16 
Cec03__ Lennard-Joes Minimum Energy Cluster [—4.4] 18 


Table 4. Results of the CEC 2019 benchmark function 


ID CDO DO GOA WOA 
Cec01 Best 53455.03 496325 170.4257 2801884356.8928 4114376128.4504 
Worst 69042.86 3532343 1917.2154 46367185532.8617 34483898 1298.6600 
Mean 103213.64 145361857875.7400 163343949975.5400 2016913382426.1000 
Std Deviation 12794.30 34144248607.8775 50100680407.7947 440383992017.2400 
Cec02 Best 18.56 18.3468 23.2134 18.3932 
Worst 18.80 18.4452 57.6369 19.1504 
Mean 19.11 18.7158 300.7714 21.2557 
Std Deviation 0.14 0.1080 60.5416 0.8178 
Cec03 Best 13.70 13.7024 13.7024 13.7024 
Worst 13.70 13.7026 13.7029 13.7025, 
Mean 13.71 13.7049 13.7062 13.7038 
Std Deviation 0.00 0.0007 0.0011 0.0003 
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Figure 3. The convergence curve of benchmark function (a) F1, (b) F2, (c) F3, (d) F4, (e) F8, (f) F9, 
(g) F10, (h) F11, () F18, (j), F19, (k) F20, (1) F21, (m)cec01, (n)) cecO2, (0)) cec03 
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Table 5. Transient response 
Algorithm Rotor Angle Output Speed Output 
Overshoot Undershoot __ Settling Time(s) _ Overshoot _ Undershoot _ Settling Time(s) 
PSS-CDO No Overshoot -0.27 785 0.0029 -0.0837 609 
PSS-DO 0.01372 -0.527 649 0.0139 -0.11 977 
PSS-GOA 0.308 -0.55 980 0.0141 -0.1287 681 
PSS-WOA 0.011 -0.7395 607 0.0257 -0.1359 543 
PSS-Conv 0.0887 -1.1244 593 0.0823 -0.1647 602 
=e i i Ty 2 If > PSS-Conv | | 
© os aj - PSS-WOA * 08 d aise 
9) Suess GOA ia —Pss-cDo 
0.15) 3f mea “| ; =--PSs-DO | 
0 100 200 300 400 500 600 700 0 500 1000 1500 2000 2500 
Time(s) Time (s) 
(a) (b) 
Figure 4. Transient response; (a) Speed response, and (b) Rotor angle response 
CONCLUSION 


This article proposes a modification of the Dandelion Optimizer method by integrating the chaotic 


algorithm and modification of the original Dandelion Optimizer (DO) to obtain optimal parameters from the 
power system stabilizer (PSS). DO is a method adopted from the movement of dandelion seeds. Dandelion is 
one of the plants that rely on wind for seed propagation. The experimental results show that the performance 
of the CDO method can increase the ability of PSS with a fully loaded system condition. PSS optimized with 
CDO can reduce overshoot speed by 96.48% and undershoot rotor angle by 75.99% compared to conventional 
methods. 
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